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Abstract 

We describe the Monte Carlo event generator WOPPER for four fermion 
production through V^-pairs including resummed leading logarithmic QED 
radiative corrections. 



Program Summary: 

• Title of program: WOPPER, Version 1.3 (April 1995) 

• Program obtainable from: crunch, ikp.physik.th-darmstadt .de in 

the directory /pub/ohl/wopper (using anonymous Internet ftp) 

• Licensing provisions: none 

• Programming language used: FORTRAN-77 

• Computer/Operating System: Any with a FORTRAN-77 environment 

'Supported by Bundesministerium fur Forschung und Technologie, Germany 
' Email: anlauf (Scrunch, ikp.physik.th-darmstadt . de 
* Email: Thorsten . OhlOPhysik . TH-Darmstadt . de 
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• Memory required to execute with typical data: 80k words 

• No. of bits in a word: 32 

• Number of program lines in distributed program, including test 
data, etc.: w 7000 (Including comments) 

• Keywords: radiative corrections, W pair production, W decays, multi- 
photon radiation 

• Nature of physical problem: Higher order leading logarithmic QED 
radiative corrections to W pair production and decay at high energy e+e~ 
colliders, including finite width of the W's. 

• Method of solution: Monte Carlo event generation 

• Restrictions on the complexity of the problem: The matrix ele- 
ments for the hard subprocess contain only the resonant contributions to 
e + e~ — > (W + W~) — > 4 fermions in the Born approximation. The final 
state fermions from the W decays are generated at the parton level. 

• Typical running time: ss 4 ms/event, depending on the energy, cuts and 
CPU. The test run takes approximately 50 CPU seconds on a 486DX2/66 
running Linux 1.2 and 35 CPU seconds on an IBM RS/6000-520 running 
AIX 3.2.5. 
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1 Introduction 



The spectacular success of the high energy electron positron colliders LEP at 
CERN and SLC at SLAC has confirmed the predictions of the Standard Model 
(SM) for the interactions between the gauge bosons and the fermions even at the 
level of electroweak radiative corrections. The high precision of the experimental 
data allows to put limits on parameters of yet unobserved particles like the Higgs 
boson through the appearance of these particles in weak loop corrections. On the 
other hand, the non-Abelian structure of the gauge sector of the SM predicts 
couplings between the electroweak gauge bosons which have not been tested 
directly. The forthcoming upgrade of LEP to LEP-II ||, [| with a center of mass 
energy of up to y/s « 200 GeV and a future 0.5-1.5 TeV e + e~ linear collider j| 
(NLC, for short) will in particular make the trilinear WWj and WW Z couplings 
directly observable through their contribution to the production of W pairs. 

The main task of LEP-II will be a precise determination of the mass and 
the width of the W and the production cross section in the threshold region |j| . 
At NLC, a measurement of the total cross section, the W angular distributions 
and in particular of their longitudinal polarization component will give a handle 
on possible anomalous couplings among the electroweak bosons induced by new 
physics beyond the SM ^, |?J , and in turn reveal an insight into the mechanism 
of electroweak symmetry breaking. 

However, possible new physics is already severely constrained by present 
LEP and SLC data, and the effects to be expected at NLC (even more so at 
LEP-II) are small. In order to extract these small effects, one has to have a 
precise knowledge of the radiative corrections within the SM. 

The electroweak radiative corrections to the production of on-shell W's to 
one-loop order are by now well established The influence of the finite 

width of the W's has been investigated in [|l0|, [Tl||. Also, the higher order 
QED corrections have been calculated in the leading log approximation (LLA) 
and their importance has been emphasized in ref. [12] . Recently an exhaustive 
overview of the standard model predictions has been provided in fl3| . 

Unfortunately, the experimental reconstruction of the W's and the deter- 
mination of their polarization is complicated by the fact that they may decay 
either into leptons with an escaping neutrino, or into hadrons, where the jet 
energies may be poorly known due to undetected particles. In addition, the 
radiative corrections due to emission of photons produce a systematic shift of 
the effective center of mass energy towards smaller values. Such effects may 
best be studied with the help of a Monte Carlo event generator. 

Although quite a few semianalytical calculations of the corrections men- 
tioned above have been available for some time, version 1.0 of WOPPER had been 
the first publicly available and complete Monte Carlo event generator. Ref- 
erence [|l4| reported a Monte Carlo generator implementation of the complete 
0(a) corrections to the production of on-shell W's. The matrix elements for 
single photon bremsstrahlung in the process e + e~ — ► (Vy+VF - ) — > 4 fermions 
and their implementation in a Monte Carlo generator is described in |T(i[ . 
but the virtual corrections are not yet included there. In the meantime, two new 
programs have been released or are about to be released in the immediate 
future 0. 

This paper describes version 1.3 of the Monte Carlo event generator WOPPER. 
It is based on the lowest order cross section for the process e + e~ — > (W + W~) — > 
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4f and focusses on QED radiative corrections in the LLA resummed to all orders 
in a and the effects of finite width of the W's. The four- moment a of the exclusive 
hard photons are generated explicitly and treated with full kinematics. The W 
decays into fermions are implemented at the parton level, including leading 
QCD corrections to the W branching fractions. 

This write-up is organized as follows: In section 2T we outline the physics 
underlying the algorithms implemented in WOPPER. The actual implementation 
is described in section ||. The parameters controlling the execution of WDPPER 
are discussed in detail in section ^ and the FORTRAN-77 interface is presented in 
section |[ Distribution notes, a listing of all external symbols and the output of 
a test run can be found in the appendices. 



2 Theoretical Background and General Features 

2.1 QED corrections at very high energies 

In the structure function formalism [jl9| |2lJ the expression for the radiatively 
corrected cross section reads 
1 



(i) 



where a is the Born level cross section of the hard process, D(x, Q 2 ) are the 
structure functions for initial state radiation, and Q 2 is the factorization scale. 

The structure functions D sum the numerically most important leading log- 
arithmic contributions 



log 



(at LEP-II and NLC energies) 



(2) 



to the electromagnetic radiative corrections to all orders. They satisfy the QED 



evolution equation |19 
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with initial condition 

D(x,m 2 e ) = 5(1 -x) . 

The solution to eq. (^|) automatically includes the very important exponentiation 
of the soft photon contributions to the radiative corrections. 

An explicitly regularized version of (^) used in the Monte Carlo implemen- 
tation is given by 

1-e 

' 9 ;D(x,Q 2 ) = |- y ^-P ee (z)D (|,<9 2 ) (6) 



a 
2^ 



dzP ee (z) 



D(x,Q 2 ). 
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It is crucial to note that the characteristics of the generated event sample do 
not depend on e, if it is chosen well below the experimental threshold for the 
detection of soft photons. Theoretically, the cross sections will remain positive 
for all values of e, but limited storage for ultra soft photons and limited floating 
point range impose a lower limit. 

The radiatively corrected cross section (|J) is implemented in a Monte Carlo 
event generator by solving the integro-differential equation (||) by iteration and 
taking into account the energy loss in the hard cross section. It is very similar 
to algorithms for quark fragmentation in QCD |p2| . As a by-product of the 
branching algorithm, the four-momenta of the radiated photons are generated 
explicitly. 

The initial state branching algorithm has already been used in the Monte 
Carlo generator KRDNDS [^3) and an improved version in the generator UNIBAB 
[ p4[ , where the implementation has been described in detail. 

The momenta of the electron and positron after initial state branching are 
then used as the input momenta for the subgenerator of the hard subprocess 
described below. 



2.2 Born Cross Section 

In the general case, there are many Feynman diagrams contributing to the 
process e + e _ — > 4 fermions at high energies, even if one requires that the quan- 
tum numbers of the final state fermions be consistent with W pair production. 
However, the contribution of the individual Feynman diagrams may be easily 
estimated by counting the number of 'resonant' propagators, where an interme- 
diate vector boson may come close to its mass shell. Using this naive estimate, 
one finds that the contribution of these 'background diagrams' is suppressed by 
a factor Tw/M\v ^2.5% for each non-resonant boson propagator, and may be 
reduced further by appropriate cuts on the invariant masses of the reconstructed 
W's. 

In the current version 1.3 of WOPPER only the Feynman diagrams with two 
resonating W's contributing to the process e~e+ — > (W + W~) — > /1/2/3/4 are 
taken into account. In fact, a full calculation finds that this approach is 
numerically justified for LEP-II energies, unless electrons are in the final state 
and no invariant mass cuts are applied. 

The implementation of the Born cross section in WOPPER is split into two 
steps, namely the production of the virtual (off-shell) W's and their subsequent 
decay into fermions including the correlations from the polarization of the in- 
termediate W's. 

This factorization in the resonant approximation can easily be seen as fol- 
lows. Since the top quark is much heavier than the W, the virtual W's will 
decay only into (almost) massless fermions. We choose the Landau gauge. The 
couplings of the would-be Goldstone bosons to the fermions are proportional to 
the fermion masses and may therefore be safely neglected. The numerator of 
the W propagators is decomposed into a sum over three physical polarizations. 

-2^ + ^ = I>m(MK(M), MM) = 0, e(k,\)-e*(k,\') = -S xx , 

(7) 
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Figure 1: The Feynman diagrams contributing to the process e + e — > W + W . 

As in ref. R], we write the amplitude for 

e-(pi)e + (p 2 ) - ^-(fc_)Fy+(fc+) - h(qi)h{q2)h{qz)h(qi) (8) 
as a coherent sum of a production and a decay amplitude for virtual Ws, 
M = M e + e -^ w+w -(pi,p 2 ;k-,\-;k + ,\ + ) (9) 



A+,A_ 



xZ(k 2 _) 1 Z{k\) 1 M w -^ flf - 2 (\^; qi ,q 2 ) M w+ _ hf - i (X + ;q 3 ,q 4: ) 



where 



z(k 2 y x 



i 



(10) 



k 2 -M^ + -£ w (k 2 ) ' 

Although we are working in the Born approximation, we shall use the renor- 
malized W self energy Y,w(k 2 ) in this expression in order to obtain a realistic 
behavior of the amplitude when a W in the intermediate state comes close to 
its mass shell. For LEP-II energies and in the on-shell renormalization scheme, 
it is sufficient to take only the leading linear term of the energy dependence of 
the imaginary part of Hw into account and to neglect the real part. 



Im E w (fc 2 ) = ^-T w {k 2 ) 



k 2 



M W T 



w± w , 



Re T, w (k 2 ) « (11) 



Taking the modulus squared of (|^) and integrating over the decay angles of 
the virtual W's, one obtains the well known resonance formula |ftoj M || for 
the total cross section for production of virtual W pairs (see figure]!). 



= / ds+ds- 



s+r w (s + ) y/sir w (s-] 



ttD(s + ) tt£)(s_) 



ct //(s;s+,s_) 



(12) 



Here a ff(s; s+, s_) denotes the off-shell cross section for the production of 
two W's, and rV(s±) is the effective decay width of the virtual VT's into light 
fermions according to eq. ([ll]). The resonance factors are essentially Breit- 
Wigner functions, but with energy-dependent width 



D(s±) ( s± - M 2 ,) 2 + S± T 2 W ( S± ) 



(13) 
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For the actual implementation it is useful to apply the mappings (7 = Fw/Myy) 

/(l + 7 2 )s± - M* \ , . 



to eq. ( 12 ) in order to get a smooth integrand suitable for a Monte Carlo rejection 
algorithm: 



1 s+ 1 s_ 



riTT -772- CT o// (s;s+,a_) 



s±=M2 .(l+ 7 tani ± )/(l+ 7 2) 



= I a( s; £+,£_) (15) 

The angular distribution of the virtual PF's in the center of mass frame is 
generated from the differential cross section da(s;s+,S—;9)/dft by a standard 
mapping and rejection algorithm. 

2.3 Coulomb singularity 

Another class of universal corrections, which is important near threshold, is the 
so-called Coulomb singularity. We implement the correction p5[ 

ff ~ = I, 1 " n arCtan (, 2/3Im/3 M ) ) {W) 

to the off shell total cross section. 



2.4 Decays of the virtual Ws 

After the momenta of the virtual Ws have been fixed, it is straightforward to 
generate the decays of the W's in the center of mass system. 

The measure of the integrations over the W-dec&y angles is independent of 
the other kinematical variables, thus leading to the simple form of eq. (|l2|). 
Hence, the decay angle distributions for fixed virtualities and production angles 
of the MK's are proportional to the modulus squared of the amplitude (||), which 
is a sufficiently smooth function of the decay angles and therefore suitable for a 
rejection algorithm. 

The final state fermions of the decay W — > //' are simply chosen with a 
probability according to the branching fractions obtained from the tree-level 
formulas. Inclusive QCD corrections to the W decays to quarks are taken into 
account to first order in as- The branching fractions for hadronic decays into u, 
d, c, s, and b quarks are given by the corresponding CKM matrix elements, 
while decays into t quarks are assumed to be kinematically forbidden because 
of |m t — Myy\ 3> r^. Currently, the fermion masses are only taken into account 
kinematically and not in the matrix elements. 

2.5 Hadronization 

If one or both of the Ws have decayed into quarks, they can optionally be 
hadronized using either the LUND string model [^?J or the HERWIG cluster 
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model Q . In the present version it is not possible^] to study color-rearrangement 
effects in purely hadronic decays. Both W decays are handled separately. 

The hadronization model can be switched at run time with the parameter 
qcdmc, which takes the values 0, 1 and 2. These correspond to no QCD MC, 
JETSET pi and HERWIG p§ respectively. 



3 Implementation of WOPPER 1.3 

Like almost all Monte Carlo event generators, WDPPER is divided into three 
parts: initialization, event generation, and termination. These are described in 
this section. For ease of use, WDPPER comes with two application interfaces, so 
that a direct call to the lower level parts in this section will never be necessary. 
These interfaces will be explained in section 0. 



3.1 Initialization 

The initializations in WDPPER are used for computing the value of variables that 
will be used frequently during event generation. Examples are the calculation of 
electroweak couplings, the maximum of the off-shell Born cross section and other 
internal steering parameters from the input parameters. This is accomplished 
by a call to the subroutine wwinit after setting the Monte Carlo parameters. 
Since WOPPER does not yet include weak corrections, most of the initializations 
performed are quite trivial. 

Finally, a standard /hepevt/ initialization record [p9| is written, which may 
be used by the analysis program. 



3.2 Event Generation 

The routine wwgen produces an event on every call. The four momenta of all 
generated particles as well as supplemental information is written to a standard 
/hepevt/ event record, where it can be read by user supplied analyzers. See 



section 3.4 below for details on WOPPER's extensive use of /hepevt/. 

The first step is the generation of the initi al s tate radiation by the branching 
routine wwbini using the algorithm of section and of the virtual masses of the 
W's occurring in the intermediate state. According to the cross section jl5|), the 
rejection weight is calculated from the ratio of the actual off-shell cross section 
to the maximum determined in the initialization step, and the event is accepted 
with a probability according to this weight. 

After the effective center of mass energy and the W virtualities have been 
fixed, the angular distribution of the virtual W's is generated in wwgt using the 
helicity amplitudes wwxhel for off-shell W pair production. Finally, the decay 
of the intermediate W's into the final state fermions is accomplished in the 
subroutine wwgdec. 



1 Unless one makes the necessary changes to the WOPPER sources in routines wwlundO and 
wwhwigO . 
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3.3 Termination 



The cross section for the generated events is obtained in the subroutine wwclos 
from the standard formula 

, # of successful trials . . 

cr t ot{s) = max { CT ( S ;£+,£_)}■————— — 17 

«';£+,£_ total # of trials 

where a is the off-shell cross section from eq. (|l5]). s' varies between and s, 
and £± in the range allowed by the specified cuts. The corresponding statistical 
error from the Monte Carlo integration is given by 

Aa tot (s) = max {*(«';£+,£_)}■ (18) 

/ (total # of trials — # of successful trials) • # of successful trials 
y (total # of trials) 3 

The cross section and the error on the cross section are placed into /hepevt/, 
where they may be read by the user-supplied analyzer. 

3.4 Additional information in /hepevt/ 

Because WOPPER uses the standard /hepevt/ event record internally, not only 
stable particles with isthep(i) = 1 will be present. Adapting the conventions 
of the HERWIG Monte Carlo ^8|, we use the following status codes 

• 101: e~ beam (positive z-direction), 

• 102: e + beam (negative z-direction), 

• 103: center of mass system of the collider, 

• 110: e + e~ hard scattering center of mass system, 

• 111: e~ before hard scattering, 

• 112: e + before hard scattering, 

• 113: virtual W~ after hard scattering, 

• 114: virtual W + after hard scattering. 

However, these entries have no physical significance and should never be used 
in any analysis (an exception to this rule are the beam particles 101 and 102, 
which are convenient for defining the reference frame and are used e.g. by the 
analyzer hepawk ]30| | for this purpose). Only the particles with status code 1 
belong to the final state as predicted by WOPPER. 

If no hadronization Monte Carlo is active, final state quarks will be entered 
as "stable" particles with status code 1. 

4 Parameters 

The parameters controlling WOPPER version 1.3 are summarized in tables |l| and 
pi They will be described in the following subsections. 
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Variable name 


semantics 


Default value 


ebeam 


e + /e~ beam energy 


250 GeV 


epol 


longitudinal e~ beam polarization 





ppol 

r r 


longitudinal e + beam polarization 





scheme 


renormalization scheme 


1 


ahpla 


/ vtf -i-J -i-^ \ 1/1/ / 


128 


ahplaO 


l/aosrifO) 


137.0359895 


massle 


m e ± 


0.51099906 • 10~ 3 GeV 


masslw 




80.22 GeV 


gammlw 


rv 


2.08 GeV 


masslz 




91.187 GeV 


gammlz 




2.492 GeV 


gf ermi 


G F 


1.16639 • 10~ 5 GeV 


sin2w 


sin 2 9w 


0.2261 


alphas 


a Q CD(M$y) 


0.12 


ckmvus 


V us = sin 9 C 


0.2196 


ckmvcb 


v cb 


0.0400 


ckmvub 


V ub 


0.0032 



Tabic 1: Physics parameters controlling WDPPER. 
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Variable name 


semantics 


Default value 


cutmin 


minimum Vt /± virtuality 


GeV 


cutmax 


maximum W ± virtuality 


y/~S = 'ZE Beam 


nevent 


Number of events 


10000 


cc 


apply canonical cuts 


.false. 


cclvl 


level of canonical cuts 





ccver 


version of canonical cuts 


1 


bstyle 


Key for QED radiation 


1 


epsiln 


Internal infrared cutoff 


10~ 5 


coulom 


Include final state coulomb corrections 


. false . 


qcdmc 


Key for QCD Monte Carlo 





rangen 


Random number generator 


1 


rseed 


Random number seed 


54217137 


rlux 


'Luxury' of ranlux 


3 


errmax 


maximum error count 


100 


verbos 


verbosity 





run id 


run identification 




stdin 


standard input 


5 


stdout 


standard output 


6 


stderr 


standard error 


6 



Table 2: Technical parameters controlling WDPPER. 
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4.1 Electroweak Parameters 



Since the present version of WDPPER does not include an electroweak library, the 
electroweak parameters, namely the masses of the electroweak bosons (masslw, 
masslz), their widths (gammlw, gammlz) and the Weinberg angle (sin2w) are 
treated as independent parameters. They enter an effective Born cross section 
and may be set directly by the user. As a special case, setting gammlw to re- 
produces the narrow- width approximation with on-shcll Ws in the intermediate 
state. 

The bulk of the non electromagnetic radiative corrections can be incorpo- 
rated into the hard cross section by using the running QED coupling ag£fl(4M^,) 
1/128 at the W scale. This value is. however, not correct for the initial state 
radiation of on-shell photons, where ctQED(0) ~ 1/137 has to be used. The 
inverse of the former value can be changed with ahpla and that of the latter 
with ahplaO. 

The presence of initial state radiation can be toggled using the parameter 
bstyle. The supported supported values are and 1, corresponding to no QED 
radiative corrections and LLA resummed initial state QED radiative corrections. 

Since version 1.3 it is possible to choose canonical input parameters for 
benchmarking LEP2 Monte Carlos |3l) by changing the value of scheme. The 
follow values are supported: 

• scheme = 0: "free scheme", all parameters are taken from the input and 
treated as independent parameters. 



scheme = 1 (default): "Gf scheme" 



• 2* KaQEpjAM^) 

sm Uw — -= — (iyi 

yf2G F Ml 



w 



w 



• scheme = 2: "sin^^r' scheme" 



V87T V n 



2 



G = ™qed{±M w , 

Tw = -^r{ 3+ ^^) (22) 

scheme = 3: "Born scheme", tree level formulas, independent sinOw, Gp, 
OiQED^Myy) = ctQED(0) and 



(23) 



G F M^ ( Iolqcd 
W = -fa " 



• scheme = -1, -2, -3: same as the positive values, except for which 
is not calculated from the standard model expression but taken from 
gammlw instead. 
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4.2 Cuts 

In the present version 1.3 of WOPPER, only cuts in the virtualities of the inter- 
mediate W's are implemented in the event generation: 

• cutmin: minimum virtuality of the intermediate W s, 

• cutmax: maximum virtuality of the intermediate W^s. 
The cuts in virtualities have to satisfy the following conditions: 

< cutmin < cutmax < \/s = 2EBeam (24) 

A value of for cutmax will automatically be reset to the available maximum, 
namely 2E Beam - 

Since version 1.3 it is possible to apply the canonical cuts for benchmarking 
LEP2 Monte Carlos |3lJ to the event record by setting cc to .true.. These 
cuts will be reflected in the caclulated total cross section. 

4.3 Monte Carlo Parameters 

The remaining, more technical Monte Carlo parameters should be almost self 
explaining. Since our branching algorithm automatically includes soft photon 
exponentiation, the results will not depend on the value of the internal infrared 
cutoff epsiln (which is measured in units of the beam energy), provided it is 
kept well below the experimental energy resolution. However, it is not advisable 
to set it many orders of magnitude lower than the default value, because this 
may result in too high photon multiplicities that will overflow internal tables. 

A note on the random number generators available: the default value 1 of 
the parameter rangen corresponds to the standard RANMAR ||^] generator, that 
has been the generator of choice for quite some time. Recently, the quality of 
the random numbers generated by RANMAR has been questioned and unwanted 
correlations have been found, that caused large systematic errors in solid state 
physics simulations |3^| . A superior variation RANLUX has been proposed f34|| , 
which is however much slower. Setting rangen to 2 switches to RANLUX which 
can be used at the 5 "luxury levels" to 4. At the highest "luxury levels", 
WOPPER will spent more than half of the total computer time in the random 
number generator. 

However, since event generation involves a lot of decisions that effectively 
randomize the subsequences used by themselves, we do not expect that the cor- 
relations in RANMAR have any significant effect on the event samples generated 
by WOPPER. The RANLUX option has been added to WOPPER for some experimen- 
tation only. It is left in only because there is no particular reason for throwing 
it out again. 

4.4 QCD Parameters 

The parameters for JETSET should be accessed trough WOPPER's lugive com- 
mand, which is translated directly to JETSET's LUGIVE () subroutine. See the 
JETSET manual (^?J for a comprehensive description of the available parameters 
and their effects. 
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WOPPER 


HERWIG 


semantics 


Default value 


hwqcdl 


QCDLAM 


Aqcd/ GeV 


0.18 


hwrms 1 


RMASS(l) 


m d / GeV 


0.32 


hwrms 2 


RMASS (2) 


m u / GeV 


0.32 


hwrms 3 


RMASS (3) 


m s /GeV 


0.5 


hwrms 4 


RMASS (4) 


mj GeV 


1.8 


hwrms 5 


RMASS (5) 


m & /GeV 


5.2 


hwrms 6 


RMASS (6) 


m t l GeV 


100.00 


hwrms 


RMASS (13) 


m e / f -/GeV 


0.75 


hwvqcu 


VQCUT 


Quark virtuality cutoff 


0.48 


hwvgcu 


VGCUT 


Gluon virtuality cutoff 


0.10 


hwvpcu 


VPCUT 


Photon virtuality cutoff 


-1.000 


hwclma 


CLMAX 


Max. cluster mass parameter 


3.35 


hwpspl 


PSPLT 


Split cluster parameter 


1.00 


hwqdiq 


QDIQK 


Max. scale for g — > diquark 


0.00 


hwpdiq 


PDIQK 


g — > diquark rate parameter 


5.00 


hwqspa 


QSPAC 


Spacelike evolution cutoff 


2.50 


hwptrm 


PTRMS 


Intrinsic px 


0.00 



"Negative values will translate to y/s. 



Table 3: HERWIG parameters accessible from WOPPER. 

Since HERWIG does not sport the equivalent of the LUGIVEO routine, its pa- 
rameters have to be accessed through the standard WOPPER access mechanisms. 
Tables || and ^ provide a list of the available HERWIG parameters and the names 
under whith they are known to WOPPER. See the HERWIG manual |2q] for a com- 
prehensive description of the effects of these parameters. 

5 FORTRAN-77 Interface 

WOPPER version 1.3 provides two application program interfaces on different lev- 
els. The higher (preferred) level consists of the command interpreter wwdcmd 
that accepts commands in form of character* (*) strings. This driver commu- 
nicates with the analyzer hepawk j3(| by default. The lower level consists of two 
FORTRAN-77 subroutine calls: wwpsrv and wopper. 
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WOPPER 


HERWIG 


semantics 


Default value 


hwipri 


IPRINT 


Printout option 


1 


hwmaxp 


MAXPR 


Max. printouts 





hwmaxe 


MAXER 


Max. errors 


10 


hwlwev 


LWEVT 


Event output unit 





hwlrsu 


LRSUD 


Sudakov input unit 





hwlwsu 


LWSUD 


Sudakov output unit 


77 


hwsudo 


SUDORD 


Sudakov order in aj 


1 


hwnrnl 


NRN(l) 


First random seed 


17673 


hwnrn2 


NRN(2) 


Second random seed 


63565 


hwazso 


AZSOFT 


Soft gluon azimuthal correlations 


.true. 


hwazsp 


AZSPIN 


Gluon spin azimuthal correlations 


.true. 


hwblli 


B1LIM 


£>-clustcr — > 1 hadron parameter 


0.0 



Table 4: Technical HERWIG parameters accessible from WOPPER. 



5.1 Higher Level Interface 

The simple commands understood by wwdcmd are (here keywords are typeset in 
typewriter font and variables in italics; vertical bars denote alternatives) 

• initialize 

Force initialization of WOPPER and write an initialization record into the 
/hepevt/ event record, which should trigger the necessary initializations 
in the analyzer. 

• generate [n] 
* wopperappl.f 

call wwdcmd ( ' init ' ) initialize the generator 

call wwdcmd ('generate 10000') generate 10000 events 

call wwdcmd ('close') cleanup 

Figure 2: Higher level FORTRAN-77 interface 
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Generate nevent events and call hepawk to analyze them. If the optional 
parameter n is supplied, nevent is set to its value. 



• close 

Write a termination record to /hepevt/, which should trigger the neces- 
sary cleanups in the analyzer. 

• statistics 

Print performance statistics (this is usually only useful for the WDPPER 
developers, who are tuning internal parameters). 

• quit 

Terminate WDPPER without writing a termination record. 

• exit|bye 

Write a termination record and terminate WOPPER. 

• set variable ival\rval 

Set physical or internal parameters. See the tables |l| and || for a compre- 
hensive listing of all variables. For example, the command ' set ahpla 127.0' 
will set the QED fine structure constant to 1/127. 

• print variable\&ll 

Print the value of physical or internal variables. Specifying the special 
variable all causes a listing of all variables known to WDPPER. 

• debug|nodebug flag 
Toggle debugging flags. 

• testran 

Test the portability of the random number generator. We use a generator 
of the Marsaglia-Zaman variety [3^| , which should give identical results 
on almost all machines. 

• banner 

Print a string identifying this version of WOPPER. 

• echo message 

Print message on standard output. 

• lugive string 

Pass string to JETSET's LUGIVE () routine. See |27) for details. 

Unique abbreviations of the keywords are accepted, i.e. ' g 1000' generates 
1000 events. The tokens are separated by blanks. Blank lines and lines starting 
with a # are ignored and may be used for comments. For portability, only the 
first 72 characters of each line are considered. 

On UNIX systems WOPPER reads the default startup files .wopper in the 
user's home directory and the current directory, if they exist. 
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* 



wwdcmd . f 

subroutine wwdcmd (cmdlin) 
character* (*) cmdlin 



else if (cmdlin. eq. 'init') 
call wopper (0) 
call hepawk ('scan') 



initialize WDPPER 
print initialization record 



else if (cmdlin. eq. 'generate') 
do 10 n = 1, nevent 



call wopper (1) 
call hepawk ('scan') 



generate an event 
analyze the event 



10 



continue 



else if (cmdlin. eq. 'close') 
call wopper (2) 
call hepawk ('scan') 



get total cross section, etc. 

finalize analysis 



else 



end 



Figure 3: Event generation loop 



5.2 Lower Level Interface 

The subroutine wopper (code) has a single integer parameter. The parameter 
code is interpreted as follows: 

• 0: initialize the generator and write an initialization record to /hepevt/. 

• 1: generate an event and store it in /hepevt/. If WOPPER has not been 
initialized yet, the necessary initializations are performed, but no initial- 
ization record is written. 

• 2: perform final calculations and write the results to /hepevt/. 

Figure ^ displays excerpts from a simplified version of wwdcmd that make the 
correspondences between the two levels of the F0RTRAN-77 interface explicit. 

WOPPER's parameters can be accessed on the lower level by the subroutine 
wwpsrv (result, action, name, type, ival, rval, dval, lval). The pa- 
rameter is specified by its (lowercase) name in the character* (*) string name. 
The string action is either 'read' or 'write' corresponding to whether the 
parameter is to be inspected or modified. The type of the parameter ( ' int ' , 
'real', 'dble', or 'lgcl') is specified in type; it is an input parameter for 
'read' and an output parameter for 'write ' . Depending on this type the value 
is passed in ival, rval, dval, or lval, respectively. The following error codes 
will be returned in the string result: ' ' : no error, 'enoarg' : invalid action, 
'enoent': no such parameter, 'enoperm': permission denied, and 'enotype': 
invalid type. 

The protection scheme implemented with this parameter handling has been 
described in p3|l . Its main purpose is to guarantee consistency of user defined 
and computed parameters in the generation phase of the Monte Carlo. 
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6 Conclusions 



We have presented version 1.3 of the Monte Carlo event generator WOPPER for 
W pair production and decays into four fermions at high energy e + e" colliders. 
The distinguishing feature of WOPPER is the inclusion of higher order electromag- 
netic corrections including soft photon exponentiation and explicit generation 
of exclusive hard photons. In contrast to fixed order calculations which have to 
be exponentiated by hand, WOPPER handles the multiphoton effects explicitly. 

The present version does not contain weak corrections. Forthcoming versions 
of the Monte Carlo generator will include weak corrections in the framework of 
effective Born cross sections |3f| . Anomalous couplings for the 7 WW and ZWW 
vertices may also be included in a future version of the generator. 
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A Distribution 

The latest release of WOPPER is available by anonymous ftp from 

crunch . ikp . physik . th-darmstadt . de 
in the directory 

pub/ ohl/ wopper 
or on the World Wide Web at the URL 

\protect\vrule widthOpt\protect\href {http : // crunch. ikp .physik. th-darmstadt .de/monte-carlos . 

Important announcements (new versions, fatal bugs, etc.) will be made through 
the mailing list 

wopper-announce@crunch . ikp . physik . th-darmstadt . de 

Subscriptions can be obtained from 

majordomoOcrunch . ikp .physik. th-darmstadt .de 

(send a message consisting of help to majordomo for instructions on how to 
subscribe, don't send such messages to the list itself). 
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B Installation 



WOPPER is distributed in PATCHY format |3fj]. Plain FDRTRAN-77 versions can 
be made available on request. Since the FDRTRAN-77 source conforms to ANSI 
X3. 9-1978, it should run without modifications on all platforms. 

B.l UNIX Systems 

On UNIX systems, the configuration, compilation and installation can be per- 
formed automatically according to the following familiar sequence: 

$ ./configure 
$ make 
$ make test 
$ make install 

Figure |J shows the command line options of the configure script for WOPPER 
on UNIX systems. This configure script has been created by the popular GNU 
Autoconfj37| package and should work on all UNIX variants. 

B.2 Non-UNIX Systems 

For non-UNIX systems configuration and compilation has to be performed man- 
ually from the CARDS file and the cradles. 

C External Symbols: Common Blocks and Sub- 
routines 

To avoid possible name clashes with other packages, all external symbols ex- 
ported by WOPPER begin with the two letters WW, except for the routine wopper 
itself and the /hepevt/ common block. 

• Common Blocks: 

The following common blocks are used by WOPPER. 

— /hepevt/, /hepspn/: standard common blocks for passing gener- 
ated events p9| . 

— /wwpcom/ : main parameter common block, holds all physical parame- 
ters. Application programs should access this common block through 
the wwpsrv routine. 

— /wwcbrn/: internal parameters used for born cross section. 

— /wwcdec/: internal parameters used for W decays. 

— /wwcevt/: internal parameters used for event generation. 

— /wwesta/: statistics. 

— /wwctri/: storage for keyword lookup. 

• Driver Program: 
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Usage: configure [options] [host] 

Options: [defaults in brackets after descriptions] 
Configuration: 



— cache-f ile=FILE 
— help 
— no-create 
— quiet, — silent 
— version 
Directory and file names 
— prefix=PREFIX 



cache test results in FILE 

print this message 

do not create output files 

do not print 'checking...' messages 

print the version of autoconf that created configure 



-exec-pref ix=PREFIX 



-srcdir=DIR 



install architecture-independent files in PREFIX 
[/usr/local] 

install architecture-dependent files in PREFIX 
[same as prefix] 

find the sources in DIR [configure dir or . .] 
— program-pref ix=PREFIX prepend PREFIX to installed program names 
— program-suf f ix=SUFFIX append SUFFIX to installed program names 
— program-transf orm-name=PROGRAM run sed PROGRAM on installed program names 
Host type: 

configure for building on BUILD [BUILD=HOST] 
configure for HOST [guessed] 
configure for TARGET [TARGET=HOST] 



— build=BUILD 
— host=H0ST 
— target=TARGET 

Features and packages: 
— disable-FEATURE 
— enable-FEATURE [=ARG] 
— with-PACKAGE [=ARG] 
— without-PACKAGE 
— x-includes=DIR 
— x-libraries=DIR 

— enable and 
— with-g77 



do not include FEATURE (same as 
include FEATURE [ARG=yes] 
use PACKAGE [ARG=yes] 
do not use PACKAGE (same as 
X include files are in DIR 
X library files are in DIR 
with options recognized: 

use GNU Fortran 77 



-enable-FEATURE=no) 



with-P ACK AGE=no ) 



-enable-verbose-patchy display all patchy output 



-enable-internal 

-enable-not ime 

-enable-pedantic 

-with- 1 ibpath=P ATH 

-with-srcpath=PATH 

-with-hepawk 

-with-jetset 

-with-herwig=CAR 

-with-cernlib 

-enable-paper-a4 

-enable-paper-us 



do not use this! 

do not use timing functions 

no IMPLICIT NONE 

use PATH for libraries 

use PATH for source files (CARs) 

use HEPAWK 

use JETSET 

use HERWIG CAR file 

use CERNLIB 

use European (A4) paper 
use US (letter) paper 



Figure 4: Comandline options of the configure script for WOPPER on UNIX 
systems. 
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— wwdriv: sample main program, which reads commands from stan- 
dard input and feeds them into wwdcmd. 

— wwdloo: command loop, reading command from a terminal or file 
and executing them. 

— wopper: the low level entry point into the generator for application 
programs. 

— wwdcmd: WDPPER's command interpreter, the preferred entry point for 
application programs. Executes a single command. 

— wwdlxi, wwdlxd, wwdlxs: Utility routines: tokenization of input. 

— wwdsig: UNIX signal handler. 

— isatty: check if this job is run interactively. 

Parameter Management: 

These routines are used to control the parameters common block /wwpcom/. 

— wwpsrv: server handling parameter changing requests. 

— wwpini: block data supplying default values. 

— wwpprn: print parameters. 

Initialization: 

— wwinit: main entry point for initializations. 

— wwigsw: initialization of electroweak parameters. 

— wwicut: initialization of internal Monte Carlo parameters. 

— wwibmx: finds maximum of Born cross section. 

— wwibn: auxiliary function for finding maximum of on-shell cross sec- 
tion. 

— wwibns: auxiliary function for finding maximum of cross section for 
symmetric virtualities. 

— wwibna: auxiliary function for finding maximum of cross section for 
one W on-shell. 

— wwibnv: 
Final calculations: 

— wwclos: calculate total cross section, errors and close the generator. 

— wwstat: statistics. 

Hard Subprocess Generation: 

— wwgen: main entry point for hard subprocess generation. 

— wwgt: generation of angular distribution of virtual W's. 

— wwgppr: generate four-momenta of virtual W's. 

— wwgdec: generate final state fermions from W decays. 

— wwdqf 1: select quark flavors in W decay. 

Branching: 
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— wwbini: generates the initial state photon radiation. 
Cross Sections: 

— wwxhel: hclicity amplitudes for W pair production. 

— wwxtot: total off-shell cross section. 

— wwxint: auxiliary function for integration over virtual W masses. 
Accessing /hepevt/: 

— wweeni: enter identification of Monte Carlo and run. 

— wweens: write summary record 

— wweent: enter one particle 

— wwenul: enter null particle 

— wwenew: clear /hepevt/ 

Hadronization: 

— wwpart: leave partons alone. 

— wwlund: JETSET interface code 

— wwhwig: HERWIG interface code 

— hwaend: HERWIG abnormal end routine 

Random Numbers: 

— wwrgen: returns a double precision uniform deviate. 

— wwrmz: random number generator RANMAR. 

— wwrlux: random number generator RANLUX. 

— wwrtst: test the portability of the random number generator. 

— wwrtmz: test RANMAR. 

— wwrtlx: test RANLUX. 

— wwrchk: used in testing RANLUX. 

— rluxat , rluxgo, rluxin, rluxut: additional RANLUX entry points 
which are not used by WOPPER. 

Utilities: 

— wwumsg: messages and error exit. 

— wwulwr: convert input to lower case. 

— wwuboo: boost a four vector. 

— wwutim: still available CPU time for this job. 

— wwuamo: multidimensional minimization. 

— wwumin: oncdimensional minimization. 

— wwuons: Gram-Schmidt procedure. 

— wwuort: another Gram-Schmidt procedure. 

Canonical cuts: 
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— adloth: apply canonical cuts in /hepevt/. 

— adloip: inner product of vectors in /hepevt/. 

— adloan: angle to beam in /hepevt/. 

— adload: add four momenta in /hepevt/. 

• Keyword search: 

(using the dynamic tries described in |$8j). 

— wwtins: insert a new keyword. 

— wwtlup: look up a (possibly abbreviated) keyword. 

— wwtnew: insert new a node into the trie. 

— wwtlen: calculate length of keyword. 

— wwtc2a: convert keyword from character* (*) to integer (*). 
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Test Run 



WOPPER version 1.3 is distributed together with a sample command file and 
hepawk script, which are given below. To run this example, the user will need to 
link WOPPER with the CERN library, because histogramming is done by HBOOK 

s 

The file sample . wopper is read from standard input (unit stdin, which is 
initialized to 5) , and sample . hepawk is read from the file SCRIPT (i.e. under MVS 
from the file which has been allocated to the DDNAME SCRIPT and under UNIX 
from the file script or from the value of the environment variable SCRIPT). 

sample .wopper 

Here is a simple WOPPER command file, setting up parameters and generating 
10000 events. 

# sample .wopper — sample WOPPER command file 

# parameters 
set ebeam 87.5 

# run 
init 

gen 10000 

close 

quit 

sample .hepawk 

This is a small hepawk script that counts the muons from the W decays and plots 
a histogram of their energy distribution. The first generated event is dumped 
to illustrate the usage of the /hepevt/ common block. 

# sample .hepawk — sample HEPAWK analyzer for WOPPER. 

BEGIN 
■C 

printf ("\nWelcome to the WOPPER test:\n"); 

printf ( 11 ***************************\n\n" ) ; 

printf ("Monte Carlo Version: */„s\n" , REV); 
printf (" Run: */.d\n" , RUN); 

printf (" Date: °/.s\n\n", DATE); 

E_max = 100; 
N_chan = 50; 

h_muon_energy 

= bookl (0, "Muon-energy" , N_chan, 0, E_beam) ; 
nr_muons =0; # initialize counter 

dumped_an_event = ; 



2G 



{ 

if (dumped_an_event == 0) 
{ 

dump ("vs") ; # Dumping first event 

dumped_an_event++ ; 

> 

for (Qp in LEPTQNS) 
{ 

if (abs(@p:id) == _pdg_muon) 
{ 

fill (h_muon_energy , @p:p:E); 
nr_muons++ ; 

} 

} 

} 

END 
{ 

# Dump some numbers 

pr intf ( " \nRESULTS : \n" ) ; 

printf ("********\n\n") ; 

printf ("Total events: °/ t d, total cross section: °/ t g pb\n" , 

NEVENT , XSECT * le9) ; 
printf ("Number of muons : */,d\n\n" , nr_muons) ; 

plotQ ; 

printf ("\ndone.\n") ; 
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sample . output 



The following output should result from the input files above, up to small round- 
off errors and different FDRTRAN-77 default output formats. 



wwdcmd 
hepawk 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 
wwigsw 



message: Starting WOPPER, Version 1.02/99, (build 0/0000) 

message: starting HEPAWK, Version 1.6 

message : ******************************************** 

message: "GF scheme" selected: 

message: Using GFERMI and ALPHA as input, calculating SIN2W. 

message: Using derived W width. 

message: Parameters used in this run: 

message: AHPLA = 128.00000 ( = l/alpha(2 M_W) ) 

message: SIN2W = .23121 (effective mixing angle) 

message: GFERMI = .116639E-04 GeV**(-2) 

message: ALPHAS = .12000 

message: GAMM1W = 2.08468 GeV (S.M. value, used) 

message: CKMVUS = .21960 CKM Matrix (Cabibbo angle) 

message: CKMVCB = .04000 CKM Matrix 

message: CKMVUB = .00320 CKM Matrix 

message : Z-e-e couplings : 

message: g_V = -.01396 

message: g_A = -.18579 

message: W-e-nu coupling: 

message: g = .23038 

message: Z-W-W coupling: 

message: g_ZWW = .57134 

message: gamma-W-W coupling: 

message: g_gWW = .31333 

message: Using energy-dependent width for W and Z propagators. 



Welcome to the WOPPER test: 




Monte Carlo Version: v01.02 (*** 00 00:00:00 1900) 
Run: 1035996352 
Date: Apr 10 03:40:00 1995 




Dumping the event record for event # 1 
There are 13 entries in this record: 


Entry # 1 is an incoming (HERWIG convention) electron 
p: ( .8750E+02; .0000E+00, .0000E+00, .8750E+02), m: 


. 0000E+00 


Entry # 2 is an incoming (HERWIG convention) positron 
p: ( .8750E+02; .0000E+00, .0000E+00, -.8750E+02), m: 


. 0000E+00 


Entry # 3 is the CMS system (HERWIG convention) 

p: ( .1750E+03; .0000E+00, .0000E+00, .0000E+00), m: 


. 1750E+03 


Entry # 4 is a null entry. 



Entry # 5 is a null entry. 
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Entry # 6 is a null entry. 



Entry # 7 is reserved for model builders. 
Entry # 8 is reserved for model builders. 



Entry # 9 is an existing photon 

p: ( .1554E+00; .6972E-03, -.1679E-03, -.1554E+00), m: 
The mother is the positron # 5. 



. 0000E+00 



Entry # 10 is an existing electron 

p: ( .5965E+02; .5462E+02, .1940E+02, .1411E+02), m: 
The first mother is the W- #7. 
The other mother is the anti-electron-neutrino # 11. 



.5110E-03 



Entry # 11 is an existing anti-electron-neutrino 
p: ( .2912E+02; -.2804E+02, -.3085E+01, -.7203E+01), m: 
The first mother is the W- #7. 
The other mother is the electron # 10. 



. 0000E+00 



Entry # 12 is an existing up-quark 

p: ( .3692E+02; .1710E+02, -.3106E+02, . 1030E+02) , m: 

The first mother is the W+ #8. 

The other mother is the anti-down-quark # 13. 



.5000E-02 



Entry # 13 is an existing anti-down-quark 
p: ( .4915E+02; -.4368E+02, .1474E+02, -.1705E+02), m: 
The first mother is the W+ #8. 
The other mother is the up-quark # 12. 



.1000E-01 



RESULTS : 
******** 



Total events: 10000, total cross section: 12.98 pb 

Number of muons : 2194 

lMuon-energy 

HB00K ID = 1 DATE 10/04/95 

160 

156 7 X 

152 XXX 

148 XXX 

144 2 XXX5 

140 X XXXX 

136 X2XXXX72 

132 2 XXXXXXXX 

128 X XXXXXXXX 

124 XXXXXXXXXX 

120 7XXXXXXXXXX 

116 XXXXXXXXXXX 
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112 


xxxxxxxxxxx 


108 


xxxxxxxxxxx 


104 


5 xxxxxxxxxxx 


100 


X7X xxxxxxxxxxx 


96 


XXX xxxxxxxxxxx 


92 


XXX2XXXXXXXXXXX 


88 


xxxxxxxxxxxxxxx 


84 


xxxxxxxxxxxxxxx 


80 


XXXXXXXXXXXXXXX7 


76 


xxxxxxxxxxxxxxxx 


72 


7XXXXXXXXXXXXXXXX 


68 


XXXXXXXXXXXXXXXXX 


64 


XXXXXXXXXXXXXXXXX 


60 


XXXXXXXXXXXXXXXXX 


56 


XXXXXXXXXXXXXXXXX 


52 


2XXXXXXXXXXXXXXXXX 


48 


xxxxxxxxxxxxxxxxxx 


44 


xxxxxxxxxxxxxxxxxx 


40 


xxxxxxxxxxxxxxxxxx 


36 


xxxxxxxxxxxxxxxxxx 


32 


xxxxxxxxxxxxxxxxxx 


28 


XXXXXXXXXXXXXXXXXX5 


24 


7XXXXXXXXXXXXXXXXXXX 


20 


xxxxxxxxxxxxxxxxxxxx 


16 


XXXXXXXXXXXXXXXXXXXX2 


12 


xxxxxxxxxxxxxxxxxxxxx 


8 


xxxxxxxxxxxxxxxxxxxxxx 


4 


7XXXXXXXXXXXXXXXXXXXXXXX! 



CHANNELS 10 
1 

CONTENTS 100 
10 
1. 

LOW-EDGE 10 
1. 

* ENTRIES = 
.00 

* BIN WID 
.10 



1 2 3 4 5 

12345678901234567890123456789012345678901234567890 

1 1 11111111111 
247090812243555433721 
34391092999413526253963823 3 1 

111112222233333444445555566666777778888899999 
2468024680246802468024680246802468024680246802468 

2194 * ALL CHANNELS = .2194E+04 * UNDERFLOW 

.2000E+01 * MEAN VALUE = .4588E+02 * R . M . S 



done. 

wwdriv: message: bye. 

D Revision History 

Version 1.3, April 1995 

• Canonical cuts and input parameters. 

• Fixed inconsistent phase conventions, which resulted in wrong angular 
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distributions. 



Version 1.2, July 1994 

• Coulomb correction. 

• Improved JETSET and HERWIG interfaces. 

Version 1.1, February 1994 

• Hadronization, JETSET and HERWIG interfaces. 

• Minor bug fixes. 

Version 1.0, 1993 

First official release, submitted to the Computer Physics Communication Li- 
brary. 
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